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ABSTRACT 
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Two  practical  detection  schemes,  the  peak-amplitude  estimator  and 
the  double  differentiator,  are  theoretically  analyzed  and  compared. 
Also,  a  formula  is  derived  which  may  be  used  to  determine  the  ef¬ 
ficacy  of  a  false-alarm  detection  system. 
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FLUCTUATIONS  IN  THE  PRESENCE  OF  GAUSSIAN  NOISE 

The  time  at  which  a  received  signal  crosses  a  certain  level  fluctuates  in  the  presence 
of  ncise.  in  this  section  theoretical  formulas  for  the  probability  density  of  the  crossing 
time  and  its  standard  deviation  are  obtained.  It  is  assumed  that  both  the  signal  and  noise 
have  been  modified  by  the  presence  of  an  intermediate  frequency  (IF)  filter,  as  shown  in 
Fig.  1  where  s(t)  is  the  signal  and  n(f)  is  the  noise  after  demodulation.  It  can  be  shown 
that  the  effect  of  the  IF  filter  passband  on  the  signal  and  noise  is  the  same  as  the  effect 
of  a  baseband  filter  of  the  same  shape.  Hence,  in  the  following  analysis,  the  IF  filter  will 
always  be  treated  an  if  it  were  a  baseband  filter.  Fig.  2  shows  the  signal  and  the  threshold 
level  a  which  must  be  crossed  to  trigger  an  output  signal  from  the  threshold  detector.  If 
an  interval  of  time  e  is  chosen  small  enough  and  if  both  s(f)  and  n(t)  are  differentiable 
functions,  then 


s(t0  +  e)  =  s(f0)  +  s'(*o>e  ] 

and  [  (1) 

n(t0  +  e)  =  n(f0)  +  n'(t0)e,  J 

where  s'(f0)  and  n'(f0)  are  time  derivatives  at  the  point  t  =  tQ.  The  probability  P  that 
s(f)  +  i(t)  crosses  the  threshold  level  during  the  interval  e  shall  now  be  calculated. 
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RECEIVER 

IF 

DEMODULATOR 

s(t>  +  n(t) 
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Fig.  1— Pulse  position  detection  system 

For  a  fixed  value  of  n'(f0),  we  must  have  s{iQ)  +  n(t0)  <  a  and  s(  t0  +  c)  +  n(f0  +  e) 
>  a.  Using  Eq.  (1),  these  two  inequalities  become  a  —  s(f0)  -  s'(t0)e  —  n'(f0)e  <  n(tQ) 

<  a  -  s(f0).  This  inequality  can  be  satisfied  only  if  n'(t0)  >  -  s'(t0).  With  the  aid  of  ' 
these  last  two  inequalities,  it  is  seen  that 

P=  I  fnn'(a,0)dadp,  (2) 

•Jas{to)-Me-pe 

where  M  =  s'(tQ)  for  notational  convenemce  and  fnn'{a,  /?)  is  the  joint  probability  density 
function  of  the  noise  and  its  derivative.  It  is  now  assumed  that  n(t)  is  a  stationary 
Gaussian  process  with  a  probability  density  function  given  by  ’ 
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/„(«)  = 


V2?o„ 


exp 


(g  -  n)2 


(3) 


TIME 

Fig.  2— Demodulator  output  signal 


where  n  is  the  mean  (expected)  value  of  n(t) 
and  o%  is  the  mean  square  of  n(t).  Since  n'(t) 
can  be  considered  the  result  of  passing  n(t) 
through  a  linear  system,  n'(t)  is  also  a  stationary 
Gaussian  process.  For  a  linear  system  with  a 
system  function  H( «)  and  a  stationary  random 
input,  it  is  easy  to  show  that 


y  «  H(0)x, 


(4) 


where  y  is  the  mean  value  of  the  output  and  x  is  the  mean  value  of  the  input.  The  sys¬ 
tem  function  of  a  differentiator  is  H(o3)  -  ju>.  Thus  it  follows  from  Eq.  (4)  that  «'  =  0. 
and  we  have  the  following  expression  for  the  probability  density  function  of  n'(t): 


4’]' 

where  is  the  mean  square  of  n'(t). 


(5) 


It  shall  new  be  shown  that  n(f)  and  n'(t)  are  uncorrelated.  For  a  stationary  process, 
we  have  the  general  relation 


*(h)n'(t2)  =  -—R(T), 


(6) 


where  r  =  — t2  and  R(r)  is  the  autocorrelation  function.  Because  R{t)  is  an  even  func¬ 

tion  for  a  real  process  n(f),  the  derivative  in  Eq.  (6)  vanishes  at  the  origin  if  it  exists 
there.  Thus 


n(t)n'(t)  =  ---R(T) 


=  0, 


7=0 


(7) 


and  n(t)  has  been  shown  to  be  uncorrelated  with  n'(t). 


Since  n(t)  and  n\t)  may  be  considered  as  the  outputs  of  two  linear  systems  with  the 
same  input  u(t),  it  follows  that  n(i)  and  n'(t)  are  jointly  normal.  Since  these  two  processes 
are  also  uncorrelated,  they  are  independent;  thus  Eq.  (2)  becomes 


-fr. 

J-M 


r**S{tO) 

'($)  W  fn(<x)da. 

~asf,to)-M  e—ffe 


(8) 


Since  e  may  be  chosen  arbitrarily  small,  Eq.  (8)  reduces  to 

P  =  efn [o  s(t0)]  f  (if  4  &fn  •(,<})  d8. 

J-M 

The  integral  may  be  evaluated  by  substitution  of  Eq.  (5)  and  a  change  of  variables.  The 
result  is 


(9) 
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P-e/Ja-sft))  [Ma*(i^)  ^C!Ip(  “|J)  ■  <10> 

where  the  subscript  on  the  time  variables  has  been  dropped  and  the  time  dependence  of 
M  is  suppressed  for  notational  brevity.  We  have  also  defined 

n  “p(-f  <u) 

At  this  point,  it  is  tempting  to  define  the  probability  density  function  for  a  threshold- 
level  crossing  at  time  t  by  fc(t)  =  P/e  in  the  limit  of  small  e.  Such  a  definition  gives  us, 
however,  a  density  function  that  cannot  be  normalized.  This  anomaly  is  due  to  the  fact 
that  there  may  be  more  than  one  threshold-level  crossing.  In  fact,  there  may  be  more 
than  one  crossing  even  when  there  is  no  signal  present  at  all. 

Crossings  that  do  not  occur  near  the  time  at  which  the  signal  level  reaches  the 
threshold  level  constitute  a  “false  alarm”  problem.  This  problem  will  be  discussed  in  the 
next  section.  For  the  present,  it  is  assumed  that  these  distant  crossings  have  boon  suitably 
excluded.  It  is  further  assumed  that  there  exists  an  interval  t1  <  f  <  t2  in  which  there  is 
a  negligible  probability  of  more  than  one  threshold  crossing  and  unity  probability  of  ex¬ 
actly  one  crossing.  We  then  can  define  a  normalized  probability  density  by  fc(t)  =  P/e  in 
this  interval  and  fc(t)  =  0  otherwise.  From  Eq.  (10),  we  have 


f c(t)  =  fnlas(t)]\M erf l^j  +-?*=  exp(-|JU  ,h<t<t2. 


fc(t)  =  0, 


where  tx  and  must  be  such  that  normalization  is  possible;  that  is, 

rf2 

fc(f)dt  =  i. 

J*i 

it  is  immediately  seen  that  and  t2  are  not  uniquely  defined  by  Eq.  (13).  To 
remedy  this  defect,  we  first  define  the  mean  threshold  time  by 


=  |  tfc(t)dt. 


We  then  require  that  the  probability  of  threshold  crossing  before  t  be  equal  to  the  prob¬ 
ability  of  threshold  crossing  after  f;  that  is,  we  require 


ff  f'2 

fc(t)dt=  I  fc(t)  dt. 


The  values  of  t^and  f2  may  be  determined  by  solving  Eqs.  fi  simultaneously. 
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A  special  case  arises  when  n  =  0,  s(f)  is  antisymmetric  about  t,  and  t  is  the  same 
time  as  that  time  when  s(f)  would  cross  the  threshold  level  in  the  absence  of  noise;  that 
is,  s(t)  =  a.  Such  a  situation  is  shown  in  Fig  2^  In  this  case  fc(t)  is  symmetric  about  f, 
and  Eq.  (15)  requires  that  fj  =  t  -  T  and  f2  =  f  +  T.  The  parameter  T  is  determined  by 
substitution  in  Eq.  (13).  We  have 


I  fc(t)dt  =  l,  (16) 

T-t 


which  may  be  solved  to  determine  T  and,  consequently,  fj  and  t2-  These  results  will  be 
referred  to  later. 

There  remains  the  task  of  justifying  the  assumption  that  led  to  Eq.  (12).  We  as¬ 
sumed  that  there  is  negligible  probability  of  more  than  one  threshold  crossing  in  the 
interval  t1  <  t  <  t2.  If  there  is  more  than  one  threshold  crossing  from  below  threshold 
to  above  it,  then  there  must  be  at  least  one  threshold  crossing  from  above  threshold  to 
below  it.  Thus  the  probability  of  multiple  threshold  crossing  in  an  interval  will  be  negligi¬ 
ble  if  the  probability  of  crossing  from  above  the  threshold  to  below  it  is  negligible.  Let 
Pr  be  the  probability  of  this  reversed  threshold  crossing  in  an  interval  between  time  t.  and 
t  +  e.  Reasoning  as  in  Eq.  (2),  it  follows  that 

rM  e-fie 

/nn'(a,W  dad&.  (17) 


Calculating  in  a  r  manner  analogous  to  that  leading  to  Eq.  (12),  we  obtain  the  probability 
density  function  of  a  reverse  crossing  in  the  interval  t.x  <t  <t2: 


(18) 


If 


0<  j  2  fr(t)  dt « 1,  (19) 

then  reverse  threshold  crossings  and  multiple  crossings  are  of  negligible  importance.  In 
this  case,  Eq.  (12)  is  a  mathematically  consistent  and  accurate  expression  for  the  threshold¬ 
crossing  density  function  in  the  vicinity  of  the  pulse  edge.  Consequently,  we  may  define 
the  standard  deviation  oc  of  the  threshold  crossing  time  by 


The  implications  of  Eq.  (19)  shall  now  be  elaborated.  From  Eqs.  (12)  and  (13),  we 

have 


r2  a  •  /  Af2  \  r*  2  /  \f\ 

eXP\2^j  =  j  M erf [a- s(f)]  dt  (21) 
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Substitution  of  Eqs.  (18)  and  (21)  into  Eq.  (19)  yields,  after  some  manipulation, 


0 <  1  -  \  Mfn[a- s(t )]  dt « 1.  (22) 

•*1 


We  conclude  that  Eq.  (22)  must  be  satisfied  for  Eqs.  (12)  and  (20)  to  be  applicable.  It 
is  noteworthy  that  Eq.  (22)  implies  that  M  >  0  in  at  least  a  subinterval  of  the  interval 
tj  <  t  <  t2.  In  particular,  if  Ilf  is  a  constant,  it  must  be  positive. 


A  special  case  of  interest  occurs  when  the  noise  is  so  weak  that  there  is  essentially 
zero  probability  of  more  than  one  threshold  crossing  in  the  interval  t1  <t  <t2.  Since 
fr(t)  =  0,  it  follows  that 


M  erf 


(23) 


The  use  of  this  relation  in  Eq.  (12)  then  gives 
fc(t)  =  Mfn[a-s(t)], 
fc(0  =  0. 


t<tlft>t2. 


(24) 


It  should  be  remembered  that  Eq.  (12)  is  valid  for  any  stationary  Gaussian  noise 
process,  including  nonwhite  Gaussian  noise.  The  formula  does  not  apply  to  non-Gaussian 
noise.  However,  it  is  shown  in  the  appendix  that  Eq.  (24)  is  approximately  valid  for  weak 
noise,  even  if  it  is  non-Gaussian. 


FALSE  ALARMS 


If  the  noise  entering  a  pulse-detecting  system  has  an  amplitude  exceeding  the  thres¬ 
hold  value,  a  spurious  detection  may  occur.  This  phenomenon  is  called  a  false  alarm.  To 
reduce  the  probability  of  a  false  alarm,  a  detector  should  include  a  rejector  circuit  which 
blocks  all  input  pulses  which  do  not  exceed  a  fixed  level  for  a  specified  period  of  time. 

let  L  be  the  rejection  level  and  T  the  time  duration  which  must  be  exceeded  if  the 
pulse  is  t  j  be  passed  by  the  rejector.  We  now  determine  the  probability  that  a  noise 
pulse  will  have  an  amplitude  greater  than  L  for  the  time  Tx.  This  probability  may  be 
deduced  by  the  same  reasoning  that  led  to  Eq.  (8).  Thus  if  ^(Tj)  is  the  probability  of 
a  noise  pulse  of  duration  Tj ,  then 


OO  JO  **)  oo 

PN(1 1)=  f  |  fn(cx)dct+  f  fn’tf)dp  f  4(a) da.  (25) 

-b  -X 


This  relation  gives  the  probability  o f  a  noise  pulse  occurring  during  an  arbitrary  interval 
of  duration  T1 .  It  is  an  approximate  relation  which  increases  in  accuracy  as  Tj  decreases. 
For  larger  Tj ,  the  expansion  of  Eq.  (1)  is  no  longer  valid  and,  consequently,  neither  is 
Eq.  (25). 
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Using  Eqs.  (3)  and  (5)  and  assuming  n  =  0,  Eq.  (25)  reduces  to 

Pw<r,)-l-i«i(|-)-  |  /„■«!)  erf df.  (26) 

Suppose  a  pulse  of  duration  T±>  T  passes  the  rejector.  The  probability  that  a  noise 
pulse  of  duration  T-±  occurred  during  the  same  interval  and  that  no  signal  was  received 
during  this  interval  is  given  by  the  joint  probability  function  Pfjt<p.  Let  Pq\n  be  the  con¬ 
ditional  probability  that  no  signal  was  received,  given  that  a  noise  pulse  occurred.  If  Pp 
is  the  probability  that  the  passed  pulse  is  a  false  alarm,  then 

Pp  =  PN,<!>  =  PNP<p\N  =  PNP<t> »  <27) 

where  the  last  equality  follows  from  the  fact  that  the  signal  and  noise  are  statistically  inde¬ 
pendent.  From  Eq.  (25)  it  is  apparent  that  Pw(Tj )  <  pn(T)  when  T1  >  T.  Since  P$  <  1, 
it  then  follows  that  Pp  <  Pjy(T).  It  can  be  shown  that  erf(x)  >  1  -  1/2  exp  (—x2/2) 
when  x  >  0.  Using  this  inequality  and  Eq.  (26), 

Pf<H1*T° 

where  T  is  the  minimum  possible  duration  of  a  passed  pulse. 


PRACTICAL  APPLICATION 

A  situation  of  common  occurrence  is  a  signal  consisting  of  pulses  which,  in  the 
vicinity  of  an  edge,  have  the  form  shown  in  Fig.  3.  The  threshold  level  a  is  set  at  half 
the  value  of  the  signal  amplitude.  In  this  symmetrical  situation  the  mean  erasing  time  is 
the  time  at  which  the  unperturbed  signal  reaches  the  threshold  level;  that  is,  s(t)  =  a.  We 
now  seek  to  evaluate  the  probability  density  function  given  by_Eq.  (12).  As  mentioned 
previously,  when  s(f )  is  antisymmetrical  about  t,  we  have  fj  =  t  —  T  and  f2  =  t  +  T.  The 
parameter  T  is  determined  from  Eq.  (16). 

If  s(f)  is  approximated  by  Fig.  3,  it  may  be  described  analytically  by 


s(f)  =  a  +  M{t  -7), 

-  a  —  a 

s(t)  =  2a, 

_  a 

t>t  +  — ’ 

M 

W 

II 

© 

,<T~r 

where  M  is  the  slope  of  the  ramp  and  is  assumed  to  be  a  constant;  that  is,  M  =  s'(t).  It  is 
also  assumed  that  T  is  large  enough  that  s(f  +  T)  =  2a;  that  is,  it  is  assumed  that  T  >  a/M. 
We  now  use  Eq.  (3),  with  ri  =  0,  and  Eq.  (29)  in  Eq.  (12).  The  result  is 


3KS 
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fc(t)  = 


exp 


M2(t-t )2 

A/ erf  | 

f  M  \ 

I  .  °n'  rrnj 

r  ^2\1 

2°« 

\°n) 

1  Vef  p' 

(  2°n  1 

—  a  —  a 
t-  —  <t<t  +  —> 
M  M 


(.-K.o-y 
4(f)  =  0,  t  <7—  T,  t  >7  +  T. 

This  equation  may  now  be  combined  with  Eqs.  (16)  and  (20),  which  yields 

Hvhvfe  “pHf  )][erf(t) "  “(t) 

-fc.  kip(-4  fe*  (’•-s*)** 


y  (30) 


(31) 


and 


_  a 

T~  w+ — a  exP| 
M  a  < 


(aog)  ‘ 1  erf(<v) 

♦sJfc-K-sf)]  RtHl; 


(32) 


The  condition  under  which  Eq.  (31)  is  a  good  approximation  is  given  by  Eq.  (22). 
Noting  that  M  =  0  for  f  >  f  +  ( a/M )  and  for  t  <  7  -  (a/M),  Eq.  (22)  becomes 


<t) 


2  erf! - 1  « 1. 


(33) 


A  tangential  signal-to-noise  ratio  is  defined  as  8.5  dB.  For  this  value  a/an  «*  1.33  and 
2  erf(-a/o„)  ~  0.18.  Thus  the  above  inequality  is  barely  satisfied  when  a  tangential  signal 
is  present.  The  left-hand  side  of  Eq^  (23)  represents  the  probability  of  a  reverse  threshold 
crossing  in  the  interval  t  —  T<t<t  +  T.  Cer¬ 
tainly,  if  this  probability  is  18%,  the  theoretical 
development  of  the  first  section  becomes  highly 
questionable.  Therefore,  it  may  be  concluded 
that  the  signal-to-noise  ratio  must  be  somewhat 
greater  than  tangential  for  Eq.  (31)  to  provide 
an  accurate  value  of  the  standard  deviation. 


To  use  Eq.  (31),  we  must  know  the  values 
of  M ,  an ,  and  on*.  Each  of  these  quantities  is  a 
function  of  the  IF  filter  shown  in  Fig.  1.  In 
addition  the  detection  process  employed  by  the 


Fig.  3— Edge  of  an  idealized  pulse 
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threshold  detector  may  alter  one  or  more  of  these  quantities.  In  the  next  two  sections, 
we  shall  consider  specific  choices  of  the  threshold  detector  and  the  IF  filter. 


PEAK-AMPLITUDE  ESTIMATOR 

A  simplified  schematic  diagram  of  a  peak  amplitude  estimator  is  shown  in  Fig.  4. 
This  circuit  is  an  adaptive  threshold  detector;  that  is,  it  will  locate  the  leading  edge  of  a 
pulse  regardless  of  its  amplitude.  In  effect,  the  threshold  level  is  adjusted  when  the  pulse 
amplitude  changes.  The  input  signal  and  noise,  s,(f)  +  n,(£),  are  fed  into  this  threshold 
detector  from  the  demodulator,  as  shown  in  Fig.  1.  The  operation  of  this  circuit  is  most 
easily  explained  if  it  is  initially  assumed  that  n^i)  =  0,  that  the  filter  bandwidth  is  wide, 
and  that  the  amplifier  has  a  gain  G  =  -2.  Suppose  that  $,•  (f)  has  the  form  indicated  in 
Fig.  5.  Referring  to  Figs.  4  and  5,  it  is  seen  that  the  output  signal  s0(f)  is  a  sharp  pulse 
rising  at  time  t  +  t,  where  t  is  the  time  delay  of  the  delay  line.  In  Fig.  5,  5  is  the  delay 
of  the  filter.  It  should  be  noticed  that  s2(t)  goes  from  above  the  threshold  to  below  it, 
whereas  our  derivations  have  assumed  the  opposite.  However,  it  may  be  easily  verified 
that  a  threshold  crossing  from  above  to  below  does  not  change  any  of  our  formulas  if  we 
merely  interpret  M  as  the  absolute  value  of  the  slope  and  a  as  the  absolute  value  of  the 
threshold  level. 


Fig.  4— Peak-amplitude  estimator 


When  noise  is  present  at  the  detector  input,  there  will  be  two  contributions  to  the 
zero-crossing  fluctuation  of  s(f).  One  contribution  is  the  noise  at  the  leading  edge  of 
s2(t);  the  other  is  the  amplitude  fluctuation  of  the  flat  portion  of  Sj  (f).  The  latter  fluctua¬ 
tion  is  due  to  the  noise  in  the  upper  branch  of  the  circuit  of  Fig.  4. 

A  filter  is  inserted  to  reduce  the  noise  in  the  upper  branch.  However,  decreasing  the 
filter  ^bandwidth  lengthens  the  rise  time  of  s1(t).  Consequently,  the  amplitude  of  Sj(t)  at 
time  t  +  r  is  decreased.  As  seen  m  Fig.  5,  this  decrease  will  cause  an  erroneous  early 
detection  of  the  pulse.  To  remedy  this  situation,  the  amplifier  of  the  lower  branch  must 
have  its  gain  adjusted.  If  the  amplifier  gain  is  G, 

s(t  +  t)  =  s1(t  +  t)  +  s2(t  +  t)  =  Sj  (f  +  t)  +  G  sfi f).  (34) 

To  ensure  that  s(f  +  r )  =  0  in  the  absence  of  noise,  it  is  necessary  that 


s  i  (t  +  t) 
i0)  ' 


(35) 
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If  Eq.  (35)  is  not  satisfied,  the  zero  crossing  of  s(f)  will  occur  at  a  time  t  +  t  +  At. 
If  At  is  small, 


Af  ~  s^_+  T)~  si(J~t'r)  +  Gsi(£) 

s'(t  +  t )  si  (F  +  r)  +  Gs.(t)  ’ 


(36) 


where  the  prime  denotes  differentiation.  This  change  in  zero-crossing  time  is  equivalent 
to  a  change  Ac  =  sj(t)Af  in  the  threshold  level  which  defines  the  pulse-edge  location. 
From  now  on,  we  assume  that  Eq.  (35)  is  satisfied  and,  consequently,  that  At  =  0. 


Since  we  are  interested  in  determining  the  position  of  the  edge  of  a  pulse,  we  can 
approximate  the  pulse  by  a  step  function  to  simplify  the  mathematics.  As  mentioned  pre¬ 
viously,  the  signal  entering  the  threshold  detector  may  be  considered  the  result  of  a  pulse 
passing  through  an  equivalent  baseband  filter  (Fig.  1).  Let  u(t)  be  the  signal  entering  this 
filter  and  h(t)  be  the  impulse  response.  Then 


«,(*)  = 


h(t-t')u(t')dt‘. 


(37) 


where  s,(f)  is  the  input  signal  in  Fig.  4.  Suppose  the  filter  transfer  function  is  H(io)  = 
\H(u))\e~lku;  that  is,  the  filter  phase  response  is  assumed  to  be  linear  over  the  range  of 
significant  values  of  |#(oo)j.  It  is  also  assumed  that  H( 0)  =  1.  If  u(t)  is  a  step  function  of 
amplitude  2a j,  occurring  at  t  -  0,  the  convolution  theorem  may  be  used  with  Eq.  (37)  to 
show  that  s,(fe)  =  a*.  Thus  we  may  make  the  identification  t  -  k.  It  then  follows  from 
Eq.  (37)  and  the  convolution  theorem  that 


L 


|H(w)|  doi. 


(38) 


\ 


It  is  immediately  noticed  that  sj(f)  is  maximum  at  t  =  t.  Referring  to  Figs.  4  and  5,  the 
absolute  value  of  the  slope  of  s(t)  at  t  =  t  +  r  is 


Af  =  |s'(f  +  t)|  =  ~s'2(t  +  t)  -  s'j  (f  +  r)  =  — Gs'-(t)  -  Sj  (t  +  r),  (39) 

where  G  is  assumed  to  be  negative.  It  is  assumed  that  the  delay  time  r  is  chosen  large 
enough  that  $j(f  +  t)  ~  0.  Then  Eqs.  (3S)  and  (39)  yield 


M  =  -G^  £^|ff(ou)|dco.  (40) 

The  mean-square  noise  associated  with  s(t)  is 


eg  =  n2(f)  =  [nj(t)  +  n2(f)]  2  =  nf(t)  +  n\{t)  +  2n1(t)n2(f).  (41) 

To  evaluate  the  right-hand  side  of  Eq.  (41),  we  let  S0(co)  denote  the  power  spectrum  of 
the  noise  entering  the  IF  filter.  Then  the  power  spectrum  of  the  noise  entering  the  thres¬ 
hold  detector  is 


Si(a>)  =  S0(co)|f7(«)l2, 


(42) 
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where  H{ co)  is  the  transfer  function  of  the  IF  filter.  The  delay  line  has  no  effect  on  the 
power  spectrum,  so 

=SS-Sf  Si(u)da>~  f  S0(m)iff(«)|2dw.  (43) 

If  we  let  Hi(co)  represent  the  transfer  function  of  the  filter  in  the  upper  branch  of  Fig.  4, 
we  similarly  obtain 

=  j  S0(o3M(a)l*lff1(u)ls  du.  (44) 

After  some  manipulation  the  last  term  in  Eq.  (41)  becomes 

oo 

^(Onjv)  =  ~  J^S0(<o)lH(uj)|2|i?1(cj)!e^r>d<<J,  (45) 

where  rr  is  the  relative  delay  time  between  the  two  branches.  In  general,  it  is  exceedingly 
difficult  to  evaluate  the  integral  of  Eq.  (45/.  However,  if  both  //(to)  and  Hx( a>)  are 
Gaussian  in  shape  and  if  Sp(o)  has  a  flat  spectrum,  then  it  is  found  that  is 

negligibly  small  if  rr  >  iT%  where  b  is  the  equivalent  “bandwidth”  of  |H(to)i2!i/1(co)|. 
Henceforth,  it  will  be  assumed  that  n1(t)n2(t)  is  negligible.  We  are  left  with 

PO  -  oo 

=  loo  SoC<^)I^(oo>12  dco  +  ~  j  50(w)lH(o:)i2|H1(to)i2  di».  (46) 

The  differentiation  of  a  noise  process  is  equivalent  to  passing  it  through  a  filter  with 
transfer  function  H0( to)  =  jco.  Using  this  fact  and  calculating  as  done  previously,  we  ob¬ 
tain  for  the  mean  square  of  n'(t) 

..  oo  cO 

o2-  =  —  f  S0(m)|//(w)|2«2dto  +  _  f  S0(co)iH(co)|2lH1  (00)120,2  (47) 

OO  •'—OO 

Once  the  filter  functions  |//(gj)|  and  |H1(gj)|  are  specified,  Sqs.  (40),  (46),  and  (47) 
may  be  used  in  Eq.  (31)  to  obtain  the  standard  deviation  of  the  threshold  crossing  time. 
As  a  specific  example,  let 


|//(w)|  =  exp 


(_J^\ 

\  2tJo  / 


If  the  bandwidth  is  defined  as  the  frequency  difference  between  the  half-power  points  of 
|//(w)|,  it  is  easy  to  show  that  the  bandwidth  B  and  the  parameter  co0  are  related  by 


\/Jog  2 ' 


(49) 


&?2rS££&3*£?s*f\ 


SfJ^SsgSiij'SjsfsSSfe! 


jsaag^^e^rtg»B?  ? 
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We  assume  that  fi1(«)  is  an  RC  “integrator”  filter;  that  is, 


1  icoRC  +  1 


(50) 


It  is  assumed  in  this  example  that  the  noise  is  white;  feat  is,  S0(a>)  -  N,  where  N  is  a 
constant.  Using  Eqs.  (48)  and  (50)  in  Eqs.  (40),  (46),  and  (47),  there  results 


M 


(51) 


L  >/2  \ 

|m 

— 

to 

- I 

A  ->0 kc]  exp 

\  ojqRC  j 

The  substitution  of  Eq.  (51)  and  the  identity  a  -  ~Gax  in  Eq.  (31)  yields  the  standard 
deviation  of  the  pulse-edge  fluctuation. 


I 

1 

l 

\ 


So  far,  we  have  ignored  the  zero-crossing-time  fluctuation  caused  by  the  zero-crossing 
detector  itself.  Usually,  this  detector  is  implemented  by  a  Schmitt-Trigger  circuit.  Theo- 
letically  the  Schmitt  Trigger  will  change  state  when  the  input  voltage  crosses  the  0-volt 
level.  In  reality,  triggering  may  occur  anywhere  within  a  voltage  range  of  2AV  volts, 
where  AV  is  tire  uncertainty  in  the  triggering  level.  If  the  absolute  value  of  the  slope  of 
the  detected  pulse  is  M  and  if  f0  is  the  mean  zero-crossing  time,  then  the  detected  zero 
crossing  may  occur  anywhere  in  the  range  f0  —  (A V/M)  <Kl0  r  (A V/M).  Defining 
c|y  as  the  mean-square  deviation  from  the  mean  zero-crossing  time  and  assuming  that 
there  is  a  uniform  probability  of  a  zero  crossing  anywhere  in  the  above  range. 


o|T 


i  (avv 

3  \  M  )  ‘ 


(52) 


This  equation  gives  the  lower  limit  of  the  standard  deviation  of  the  time  fluctuation  of  a 
pulse  edge;  that  is,  this  equation  gives  the  standard  deviation  in  the  total  absence  of  noise 
and  under  ideal  operating  conditions  elsewhere  in  the  detector. 


DOUBLE  DIFFERENTIATOR 

It  is  observed  from  Eq.  (46)  that  the  presence  of  the  upper  branch  of  Fig.  4  causes 
an  increase  in  the  noise  power  at  the  zero-crossing  detector.  It  is  worthwhile,  therefore, 
to  investigate  a  detection  system  which  contains  only  a  single  branch.  Such  a  system  is 
illustrated  in  Fig.  6,  where  s,(t)  represents  the  pulse  after  passage  through  the  IF  filter 
and  s,-(f)  is  determined  by  Eq.  (37).  Once  again,  the  input  to  the  IF  filter  is  assumed  te- 
be  a  step  function  of  amplitude  2a;  the  filter  transfer  function  is  H( co)  =  |f/(<o)!e~-/’w*.  It 
follows  from  Eq.  (38)  that  after  passage  through  the  double  differentiator,  the  signal  is 


s(0  =  sl"(t)=^  I  ;'wlH(w)ierw(*-0  doi. 


(53) 


i 

i 

*£ 

JA 
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sj(t)  +  n(U) 

~  1  —  . . 

DOUBLE 

1 

»(»>  +  r.(t> 

ZERO-CROSSIto 

1 

OUTPUT 

DIFFEREh  HATOR 

DETECTOR 

Fig.  6— Denble  differentia  toy  system 


It  Is  seen  that  s(t)  =  0  at  t  =  t,  which  is  the  half-amplitude  point  of  s3(f)-  Thus  the 
double  differentiator  has  coaverted  the  half-amplitude  time  to  the  zero-crossing  time. 
When  Gaussian  noise  is  present,  we  may  apply  Eq.  (31)  once  the  values  of  Af,  o|,  and 
ot2«  are  known  for  tte  system  of  Kg.  6.  From  Eq.  (53) 

CO 

Af  -  l$'(t)I  =  —  f  w2iff(o>)|  dw  (54) 

®  J—oo 


iSrce  the  power  spectrum  of  n(:)  js  given  by  S(c*>)  =  S0(cj)|ff(co)i2to4,  it  follows  that 


and 


= ^  s°(wi^w*'2w4  *°» 
3"'  =  i  f  S0(«)|//(W)|2W6  dtO, 

v—on 


(55) 


rthere  o0(a>)  u  the  power  spectrum  of  the  noise  entering  the  I?  ttter.  If  the  noise  is 
white  Gaussian,  then  S0(<*;/  =  N,  a  constant.  If  the  filter  transfer  function  is  Gaussian. 
then  Eq.  (48)  may  be  applied  to  Eqs.  (54)  and  (55)  with  the  result  that 


and 


J 


(56) 


In  the  general  case,  a  comparison  between  the  double  differentiator  and  the  peak- 
amplitude  estimator  is  difficult  without  the  aid  of  a  computer.  However,  in  the  sper  ial 
case  of  weak  noise,  we  may  use  Eq.  (24)  to  obtain  an  approximate  comparison.  Wc  as¬ 
sume  that  the  signal  s(f)  has  the  form  shown  in  Fig.  3  and  is  described  by  Eq.  (28).  We 
now  combine  Eqs.  (24),  (29),  (3),  and  (20).  Since  n  ~  d  and  M  -  0  outside  the  interval 
t  ~  (c/M)  (a/Af),  we  calculate 


(57) 


feivXito;  tauAubfc  fAiw  >!!«(  v  •* 
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The  threshold  value  is,  as  usual,  equal  to  ->rte  half  the  amplitude  of  the  signal.  When  the 
iitMse  is  weak,  the  signal-to-noise  ratio  is  high.  Thus  a/on  has  a  large  value,  and  the 
hracketed  term  in  Eq.  (57)  has  a  value  slightly  les  than  unity.  We  conclude  that 


/t2  <  —O. 

°c 


(58) 


We  now  evaluate  Eq.  (58)  for  the  double  differentiator.  From  Eqs.  (56)  and  (58) 


( °c)dD  < 


16fl2O30  ' 


(59) 


The  corresponding  equation  for  the  peak-amplitude  estimator  is  obtained  from  Eq.  (51). 
We  have 


(o\ })pa  > 


y/nN 

4afcj0' 


(60) 


To  make  a  fair  comparison,  we  must  assume  that  both  detectors  are  fed  by  signals  of 
equal  amplitude  and  that  both  detectors  have  half-amplitude  threshold  levels.  With  this 
assumption,  a1  =  a,  and  Eqs.  (59)  and  (60)  may  be  combined.  The  result  is 


<f(°f)iM-  (61) 

Thus  the  double  differentiator  gives  a  significantly  superior  performance  with  respect  to 
the  standard  deviation  of  the  pulse-edge  location.  The  only  drawback  is  that  a  satisfactory 
double  differentiator  is  difficult  to  implement. 
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APPENDIX 


Non-Gaussian  Noise 


The  preceding  derivations  have  all  assumed  that  the  noise  is  Gaussian.  In  this  sec¬ 
tion,  the  noise  is  assumed  to  be  weak  but  not  necessarily  Gaussian.  If  the  noise  is  weak, 
it  is  intuitively  reasonable  to  assume  that  there  is  only  one  threshold-level  crossing  in  the 
vicinity  of  the  signal  edge.  Formally,  it  is  assumed  that  the  probability  of  a  reverse  cross¬ 
ing  (from  above  the  threshold  to  below  it)  is  zero. 

Let  Fn(x)  be  the  probability  distribution  function  of  the  noise;  that  is,  Fn(x)  is  equal 
to  the  probability  that  n(t)  <  x.  Once  again,  a  is  the  value  of  the  threshold  level  and  s(t) 
represents  the  signal.  Since  there  can  be  only  one  crossing,  the  probability  that  no  cross¬ 
ing  occurred  before  the  time  t  is  equal  to  Fn  [a  -  s(f)] ,  for  if  n(t)  +  s(f)  <  a,  then  no 
crossing  could  have  occurred  before  the  time  t.  Let  Fc(t)  be  the  probability  distribution 
function  of  a  crossing;  that  is,  Fc(t)  is  equal  to  the  probability  that  a  crossing  occurred 
before  the  time  f.  From  this  definition  and  the  preceding  statements, 

l-Fc(t)  =  Fn[a-s(f)].  (Al) 

Differentiation  of  a  distribution  function  gives  a  density  function.  Thus  at  all  points 
where  a  derivative  is  defined,  we  have  the  following  density  function  for  a  threshold 
crossing: 


fc(t)  =  s'(t)fn[a-s(t)  ]  (A2) 

In  general,  this  expression  pertains  to  the  vicinity  of  the  signal  edge,  where  it  has  been 
assumed  that  only  one  threshold-level  crossing  occurs.  Using  the  same  normalization  pro¬ 
cedures  as  in  the  first  section,  Eq.  (A2)  leads  to  Eq.  (24),  where  M  =  s'(t )  by  definition. 
Originally,  Eq.  (24)  was  derived  by  assuming  a  single  crossing  and  Gaussian  noise.  The 
present  derivation  shows  that  Eq.  (24)  is  valid  for  any  noise  statistics,  so  long  as  the 
angle-crossing  assumption  applies. 

It  is  noteworthy  that  Eq.  (12)  reduces  to  Eq.  (24)  when  an'  -*  0,  which  occurs  when 
the  noise  is  weak.  This  observation  is  a  manifestation  of  the  connection  between  the 
single-crossing  assumption  and  weak  noise. 


